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Abstract 

We consider super symmetric theories where the gravitino is the hghtest superpar- 
ticle (LSP). Assuming that the long-hved next-to-hghtest superparticle (NSP) is a 
charged slepton, we investigate two complementary ways to prove the existence of 
supergravity in nature. The first is based on the NSP hfetime which in supergrav- 
ity depends only on the Planck scale and the NSP and gravitino masses. With the 
gravitino mass inferred from kinematics, the measurement of the NSP lifetime will 
test an unequivocal prediction of supergravity. The second way makes use of the 
3-body NSP decay. The angular and energy distributions and the polarizations of 
the final state photon and lepton carry the information on the spin of the gravitino 
and on its couplings to matter and radiation. 



1 Introduction 



Deciphering hidden symmetries in nature has been one of the most exciting and challeng- 
ing tasks in physics. Most recently, the discovery of the massive W and Z gauge bosons 
has established a spontaneously broken gauge symmetry as the basis of the electroweak 
theory. Here, we discuss how one may discover the massive gravitino, which would es- 
tablish spontaneously broken local supersymmetry as a fundamental, hidden symmetry 
of nature. 

If the theory underlying the standard model is supersymmetric, one may find su- 
perpartners of quarks, leptons and gauge bosons at the Tevatron, the LHC or a future 
Linear Collider. Even though an exciting discovery, this would still not answer the ques- 
tion how supersymmetry is realized in nature. To identify supersymmetry as an exact, 
spontaneously broken symmetry requires evidence for the goldstino. Only the discov- 
ery of the massive spin-3/2 gravitino, containing the spin-1/2 goldstino, would establish 
supergravity 1] with local supersymmetry as the fundamental structure. 

In general it is difficult to detect gravitinos since their couplings are Planck scale 
suppressed. However, evidence for the gravitino may be obtained in collider experiments 
if it is the lightest superparticle (LSP). The gravitino mass may be of the same order as 
other superparticle masses, hke in gaugino mediation j2] or gravity mediation P|. But 
it might also be much smaller as in gauge mediation scenarios jl]. As we shall see, a 
discovery of the gravitino appears feasible for gravitino masses in the range from about 
1 GeV to 100 GeV. As LSP the gravitino is also a natural dark matter candidate. 

We will assume that the next-to-lightest superparticle (NSP) is a charged slepton. 
This is a natural possibility with respect to the renormalization group analysis of super- 
symmetry breaking parameters. Scalar leptons may be produced at the Tevatron, the 
LHC or a Linear Collider. They can be directly produced in pairs or in cascade decays of 
heavier superparticles. The NSP lifetime is generally large because of the small, Planck 
scale suppressed coupling to the gravitino LSP. 

The production of charged long-lived heavy particles at colliders is an exciting pos- 
sibility (5j. Sufficiently slow, strongly ionizing sleptons will be stopped within the de- 
tector. One may also be able to collect faster sleptons in a storage ring. In this way 
it may become possible to study NSP decays. The dominant NSP decay channel is 
£ — > £ + missing energy, where i and i denote slepton and lepton, respectively. 

In the following we shall study how to identify the gravitino as cause of the missing 
energy. First, one will measure the NSP lifetime. Since the gravitino couplings are fixed 
by symmetry, the lifetime is predicted by supergravity given the gravitino mass, which 
can be inferred from kinematics. In a second step spin and couplings of the gravitino or 
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the goldstino can be determined from an analysis of the 3-body decay £ — » £ + ■?/'3/2 + 7. 



2 Gravitino mass 

To be specific, we focus in the following on the case where the scalar lepton r, the 
superpartner of the r-lepton, is the NSP. It is straightforward to extend the discussion 
to the case where another scalar lepton is the NSP. As we shall see, phenomenologically 
particularly interesting is the case where the gravitino is not ultra-light, which implies a 
long NSP lifetime. 

At LHC one expects O{10^) NSPs per year which are mainly produced in cascade 
decays of squarks and gluinos The NSPs are mostly produced in the forward di- 
rection [ZI which should make it easier to accumulate rs in a storage ring. In a Linear 
Collider an integrated luminosity of 500 fb~^ will yield C(IO^) rs [8 . Note that, in a 
Linear Collider, one can tune the velocity of the produced rs by adjusting the e'^e~ 
center-of-mass energy. 

A detailed study of the possibilities to accumulate r NSPs is beyond the scope of 
this Letter. In the following we shall assume that a sufficiently large number of rs can 
be produced and collected. Studying their decays will yield important information on 
the nature of the LSP. In the context of models with gauge mediated supersymmetry 
breaking the production of r NSPs has previously studied for the Tevatron .9j, for the 
LHC [ini and for a Linear Collider 

The NSP r is in general a linear combination of rR and r^, the superpartners of the 
right-handed and left-handed r-leptons r^ and rL, respectively, 

r = cos(v?r) + sin(v9^)rL . (1) 

The interaction of the gravitino with scalar and fermionic r-leptons is described by 
the lagrangian [T2], 



where Di, rR = {du + ieA^)T^. Here denotes the gauge boson, and Mp = (SttGn) 
is the reduced Planck mass. The interaction lagrangian of tl has an analogous form. 

The r decay rate is dominated by the two-body decay into r and gravitino. 
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where rrir = 1.78 GeV is the r mass, rrir is the r mass, and m3/2 is the gravitino mass. 
Neglecting rrir, we arrive at 

p2-body ^ 1^4 I ^ 2l± 1 (4) 

For instance, = 150 GeV would imply a lifetime of ^ 78 s or F-^ ^ 4.4 y for a 
gravitino mass of 7/13/2 = 0.1 GeV or 171^/2 = 75 GeV, respectively. The crucial point is 
that the decay rate is completely determined by the masses and m3/2, independently 
of other SUSY parameters, gauge and Yukawa couplings. The mass nir of the NSP will be 
measured in the process of accumulation. Although the outgoing gravitino is not directly 
measurable, its mass can also be inferred kinematically unless it is too small, 

= + ~ '2.m:fEj. . (5) 

Therefore, the gravitino mass can be determined with the same accuracy as Er and m^, 
i.e. with an uncertainty of a few GeV. 

Comparing the decay rate (jH]), using the kinematically determined 7713/2, with the 
observed decay rate, it is possible to test an important supergravity prediction. In other 
words, one can determine the 'supergravity Planck scale' from the NSP decay rate which 
yields, up to 0{a) corrections. 
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The result can be compared with the Planck scale of Einstein gravity, i.e. Newton's 
constant determined by macroscopic measurements, Gn = 6.707(10) ■ 10~^^ GeV~^ 



M^(gravity) = (87r Gn)"= (2.436(2) ■ 10^^ GeV)' . (7) 

The consistency of the microscopic and macroscopic determinations of the Planck scale 
is a crucial test of supergravity. 

Furthermore, the measurement of the gravitino mass yields another important quan- 
tity in supergravity, namely the mass scale of spontaneous supersymmetry breaking. 



MsusY = V V3Mp 7773/2 . (8) 

This is the analogue of the Higgs vacuum expectation value v in the electroweak theory, 
where v = V^mw/g = {2V2Gf)-^/^. 
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3 Gravitino spin 



If the measured decay rate and the kinematically determined mass of the invisible particle 
are consistent, we already have strong evidence for supergravity and the gravitino LSP. 
In this section we analyze how to determine the second crucial observable, the spin of 
the invisible particle. 

To this end, we consider the 3-body decay tr — »■ tr + V'3/2 + 1, leaving final states 
with W- or Z-bosons for future studies. We only consider the diagrams of Fig. Q and 
restrict ourselves to a pure 'right-handed' NSP tr. Here, we have neglected diagrams with 
neutralino intermediate states, assuming that they are suppressed by a large neutralino 
mass. 



7 7 7 




■03/2 "03/2 

(a) (b) (c) 

Figure 1: Diagrams contributing to r ^ r + V'3/2 + 7 at tree level. We do not take into account 
the diagram with a neutralino intermediate state. It turns out that (a) is the crucial ingredient 
to prove the spin-3/2 nature of the gravitino. 

In order to prove the spin-3/2 nature of the invisible particle, we compare the 3-body 
decay with final state gravitino with the corresponding decay involving a hypothetical 
neutralino A. As an example, we consider the Yukawa coupling^, 

^Yukawa = k {t^\ T + T^X t) + h.C. . (9) 

Accidentally, the coupling h could be very small, such that the r decay rate would be 
consistent with the rate given in Eq. (jH)). 

Also the goldstino has Yukawa couplings of the type given in Eq. 0. The full inter- 
action lagrangian is obtained by performing the substitution ^pf^ V^m^^M^ 

^This interaction would arise from gauging the anomaly free U(l) symmetry Lr — L^^, the difference 
of r- and /i-numbcr, in the MSSM, with A being the gaugino. 
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supergravity lagrangian. Using the equations of motion one finds for the non- derivative 
form of the effective lagrangian |14j . 
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where we have neglected a quartic interaction term which is irrelevant for our discus- 
sion. Note that the goldstino coupling to the photon supermultiplet is proportional to 
the photino mass rriTy. As a consequence, the contribution to r-decay with intermediate 
photino is not suppressed for very large photino masses. As we shall see, this leads to 
significant differences between the angular distributions for pure Yukawa and goldstino 
couplings. 




(a) Kinematical configuration. The 
arrows denote the momenta. 




(b) Characteristic spin-3/2 process. The 
thick arrows represent the spins. 



Figure 2: (a) shows the kinematical configuration of the 3-body decay, (b) illustrates the char- 
acteristic spin-3/2 process: if photon and r-lepton move in opposite directions and the spins 
add up to 3/2, the invisible particle also has spin 3/2. 

In the following we discuss two methods to determine the spin of the invisible particle. 
The first one is based on a double differential angular and energy distribution, the second 
one makes use of the angular distribution of polarized photons. 

In r-decay both, photon and r-lepton will mostly be very energetic. Hence the photon 
energy Ej and the angle 6 between r and 7 can both be well measured (cf. Fig. |2(a)| ). 
We can then compare the differential decay rate 

1 d2r(r^r + 7 + X) , , 

A{E„cose) = — ^ ^ , 11 

a T:f dE^ d cos 

for the gravitino LSP (X = 1P3/2) and the hypothetical neutralino {X = A) with pure 
Yukawa coupling. Details of the calculation are given in appendix El 
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(a) Gravitino (b) Spin-1/2 neutralino 



Figure 3: Contour plots of the differential decay rates for (a) gravitino V'3/2 (b) neutralino 
A. m;^ = 150 GeV and mx = 75 GeV [X = A). The boundaries of the different gray shaded 
regions (from bottom to top) correspond to A{E-y, cos 9)[GeV'^] = 10^^, 2 x 10^^, 3 x 10^-^, 4 x 
10~^,5 X 10~^. Darker shading implies larger rate. 

In the forward direction, cos 6' > 0, bremsstrahlung (cf. Fig. |l(b)| dominates, and 
final states with gravitino and neutralino look very similar. In the backward direction, 
cos^ < 0, the direct coupling Fig. |l(a)| is important, and the angular and energy distri- 
bution differs significantly for gravitino and neutralino. This is demonstrated by Fig. El 
where = 150 GeV and 7713/2 = 75 GeV {mx = 75 GeV). The two differential distribu- 
tions are qualitatively different and should allow to distinguish experimentally gravitino 
and neutralino. 

Note that even for very small masses 777,3/2 and mx, the differential decay rates A for 
gravitino ?/'3/2 and neutralino A are distinguishable. In this small mass limit, '03/2 can 
effectively be described by the goldstino x (with the interaction (fTUIl ). and the differential 
decay rates for ip3/2 and x essentially coincide. The discrepancy between A and x stems 
from the additional photino contribution, as discussed below (fTIUl . This makes it possible 
to discriminate the goldstino from the neutralino even for very small masses. 

The second method to test the spin-3/2 nature is intuitively more straightforward 
though experimentally even more challenging than the first one. The main point is ob- 
vious from Fig. |2(b)| where a left-handed photon and a right-handed r move in opposite 
directions. Clearly, this configuration is allowed for an invisible spin-3/2 gravitino but 
forbidden for spin-1/2 neutralino. Unfortunately, measuring the polarizations is a difficult 
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task. 

As Fig. |2(b)| illustrates, the spin of the invisible particle influences the angular distri- 
bution of final states with polarized photons and r-leptons. Again the difference between 
gravitino and neutralino is due to the direct coupling shown in Fig. 1(a) and most sig- 
nificant in the backward direction. An appropriate observable is the angular asymmetry 



dr ,~ dr 

[tr ^ tr 7r a) - -(tr ^ tr 7l a) 



A (r.r^^a\ — d cos 9 dcos^ (Tn\ 

^RLlcosfc^j - — 2 dT — I ' ^ ^ 

3 -^{tr -^Tr + 'Jr + X) + -(tl Tr + 7l + X) 

d cos t) a cos u 

where X denotes gravitino (X = 1/^3/2) and neutralino (X = A). Here, we also study the 
angular asymmetry for a pseudo-goldstino in the final state {X = x)- Like a pseudo- 
Goldstone boson, the pseudo-goldstino has goldstino couplings and a mass which explic- 
itly breaks global supersymmetry. Notice that, as mentioned above, the photino does not 
decouple in this case. 

The three angular asymmetries are shown in Fig. |3]for = 150 GeV and different 
masses of the invisible particle. As expected, the decay into right-handed r and left- 
handed photon at = vr is forbidden for spin-1/2 invisible particles (x and A), whereas it 
is allowed for the spin-3/2 gravitino. This is clearly visible in Figs. ^c,d); for small grav- 
itino masses the goldstino component dominates the gravitino interaction as illustrated 
by Figs.|li:a,b). 

Our discussion is easily generalized to the case where the NSP is a linear combination 
of Tr and tl. One then needs further information on the left-right mixing angle (pr, which 
could be provided by a direct measurement of the r-polarization or by the coupling to 
VT-boson. 

Let us finally comment on the experimental feasibility of the gravitino spin deter- 
mination. The angular distribution of the 3-body decay is peaked in forward direction 
{6 = 0). Hence, a large number of events is needed for the spin measurement. Com- 
pared to the 2-body decay, backward (cos6' < 0) 3-body decays are suppressed by 
~ 10^^ X a ^ 10~^. Requiring 10. . . 100 events for a signal one therefore needs 10^ 
to 10^ rs, which appears possible at the LHC and also at a Linear Collider according to 
the discussion in Sec. EI 

4 Gravitino cosmology 

The existence of gravitinos imposes severe constraints on the early history of our universe. 
If the gravitino is the LSP and stable, as assumed in our investigation, there are two 
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Figure 4: Angular asymmetries for gravitino "03/2 (solid curve), goldstino x (dashed curve) and 
neutralino A (dotted curve). We use m;^ = 150 GeV and cut the photon energy below 10% of 
the maximal photon energy (cf. App.^. Note that the asymmetries only depend on the ratio 

r = m\/ml {X = V'3/2, X, A). 

important constraints which we shall briefly discuss in this section. The first one arises 
for large reheating temperatures Tr after inflation, which may lead to a relic gravitino 
abundance exceeding the observed cold dark matter density. This 'overclosure' constraint 
implies an upper bound on the reheating temperature Note, however, that there 
are several mechanisms which avoid this constraint and which, in addition, explain the 
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observed cold dark matter in terms of gravitinos [THj . 

The second constraint concerns the decay of the long-hved NSP r. If it occurs during 
or after nucleosynthesis (BBN), it may spoil the successful predictions of BBN |151ll7j . A 
recent detailed analysis jTH] shows that the hadronic decay of a heavy particle during or 
after BBN imposes severe constraints on its abundance and lifetime. If the 3-body decay 
T —>■ ^/'3/2 + r + Z is allowed, one finds the upper bound on the r-lifetime (F^)^^ < 10^ s, 
or equivalently, 1123/2 < 0.4 GeV (m^/lSO GeV)^/^. Note that in the case of non-zero 
left-right mixing ipr, processes involving W also have to be taken into account. On 
the other hand, if hadronic r decays are sufficiently suppressed, which is the case for 
rrir — "^3/2 < ^z, only the effect of the electromagnetic NSP decay has to be taken 
into account. The allowed mass range is then extended to 100 GeV < rrir ^130 GeV and 
{rrir — mz) ^ ^ 35 GeV for a typical pair annihilation cross section cr^ of rs, and 
to 100 GeV < < 350 GeV and (m^ — mz) < "^3/2 ^ 260 GeV, if is enhanced by a 
factor 100 [20] • Note that for larger gravitino masses the spin determination is easier, as 
discussed in Sec. El 

Finally, we should mention that the above BBN constraints disappear if there is 
sufficient entropy production between the decoupling of the NSP at Td ^ C(lOGeV) 
and BBN at Tbbn — C(MeV), or if the reheating temperature Tr is lower than T^, 
^BBN <Tr < Td, so that there is no thermal production of rs. 

5 Conclusions 

We have discussed how one may discover the massive gravitino, and thereby supergravity, 
at the LHC or a future Linear Collider, if the gravitino is the LSP and a charged slepton 
is the NSP. With the gravitino mass inferred from kinematics, the measurement of the 
NSP lifetime will test an unequivocal prediction of supergravity. For gravitino masses 
^3/2 ^ 0{GeV) also the determination of the gravitino spin appears feasible. 
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A 3-body slepton decays 

This appendix provides some details of the calculations leading to the results of Sec. El 
For simplicity we take m,- = 0. The diagrams of Fig. Q yield 

J2 M{rn ^ rn + IL + ^3/2) = -—f {I - r - 2r^) ^- ' \ \ '- ^ 

(A.la) 

2e^m\ r (l-r)2(2-z)(l-r + 2r/) 



spins 



4r zr]^ 

(1 - rf\2 - {\ - r) z\ 12 - 10 (1 - r) ^ - 2 r (4 - r) ^2 ^ (1 + r - 2 r^) 

Tz{\ — r — 2r]) rz'^{l — zri) 

[2 - {l-r)zf [-8 + 6(1 -r) 2 - (1 - 7r) z"^] -2r][2z-{l -4r)^2] 
r z"^ {1 — zrjy r z"^ 



(A.lb) 



where r = Tn^i^/rn^, rj = E^/m^ and z = 1 — cos 6*. The corresponding transition 
probabilities for the (hypothetical) spin-1/2 particle A read (now with r = m|/m|) 

V - TR + 7L + = 2a) 

spins 

,m9 , , , (1 -r)2(2 - 2) +4r2r72 

- rn + 7R + A p = e' h' ^ ^ ^ ^ . / • A.2b 

2 2; rv'^ll — r — 2 771 

spins ^ ' 

In the case of pseudo-goldstino x with interactions described by Eq. (jlUj) one has to 
include the diagram with the photino intermediate state. We then obtain in the limit of 
a large photino mass rriTy (with r = rn^/m~) 

E, . ./~ m9 9 (2 — z)(l — r — 2 r?) 
|A1(m-.K + 7L + x)P = e^3^^^ \ (A.3a) 

spins d/z r 

El . m9 9 ""^y \ il — rYi2 — z) + Ar zrf 

^2..,^(l-.-2,)|2-(l-r).] ^ 2,-(l-r-2,) _ _ (^,3,) 

(1 — z r/)"' 1 — ZT] J 

The limit r — > in Eqs. ()A.3|) yields the results for massless goldstino. Indeed, as can be 
seen by straightforward calculation, they precisely reproduce the massless limits of grav- 
itino transition probabilities, which are obtained by taking the limit r — in Eqs. ()A.1|) 
while keeping the SUSY breaking parameter -\/3'^3/2^p finite. 
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In order to present the angular distribution, we perform an integration, 



El 



max 



dr 1 f Er [rrir {rrir - 2 E^) - m^j 2 



, - , , dE^ " ^ " ^ " — ^ , (A.4) 

dcos^ 1287r^ J mr [mr- {1- cose) Erf 

where mx denotes the mass of the invisible particle. Notice that the |A^ps have singular- 
ities coming from soft photon oc l/f]), soft r {1/Ej. oc l/[l — r — 2ri]) and a coUinear 
divergence (1/[1 — cos^] = 1/z). The last two are not really divergent for finite m^. We 
remove the soft photons from the rates. This procedure is justified by the limited reso- 
lution of real detectors. The requirement that E^ > 6 E^^x (where -Emax = rrir^l — r)/2 
is the maximal energy of the photon) leads to 



nir {m~ — rn\) (1 — 5) 
In Sec. El 5 = 0.1 is used, i.e. the photons are cut below 10% of their maximal energy. 



Er < Er-{5) = , ' \ \ ^ . (A.5) 

2mi — d (mi — mir) (1 — cosb*) 
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